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ABSTRACT 

The accretion disk around a rotating magnetic star (neutron star, white dwarf or T Tauri star) is 
subjected to periodic vertical magnetic forces from the star, with the forcing frequency equal to the 
stellar spin frequency or twice the spin frequency. This gives rise bending waves in the disk that may 
influence the variabilities of the system. We study the excitation, propagation and dissipation of these 
waves using a hydrodynamical model coupled with a generic model description of the magnetic forces. 
The m = 1 bending waves are excited at the Lindblad/vertical resonance, and propagate either to larger 
radii or inward toward the corotation resonance where dissipation takes place. While the resonant torque 
is negligible compared to the accretion torque, the wave nevertheless may reach appreciable amplitude 
and can cause or modulate flux variabilities from the system. We discuss applications of our result to 
the observed quasi-periodic oscillations from various systems, in particular neutron star low-mass X-ray 
binaries. 

Subject headings: accretion, accretion disks - hydrodynamics - waves - binaries: general - stars: 
magnetic helds - stars: neutron 



1. INTRODUCTION 

Disk accretion onto magnetic central objects occurs in 
a variety of astrophysical systems, ranging from classical 
T Tauri stars (e.g., Bouvier et al. 2007a), cataclysmic 
variables (intermediate polars; e.g. Warner 2004), to 
accretion-powered X-ray pulsars (e.g. Lewin & van der 
Klis 2006). The basic picture of the disk-magnetosphere 
interaction is well known: The stellar magnetic field dis- 
rupts the accretion flow at the magnetospheric boundary 
and funnels the plasma onto the polar caps of the star 
or ejects it to infinity. The magnetosphere boundary is 
located where the magnetic and plasma stresses balance, 
r m = ^[^/{GMM 2 )) 1 / 7 , (1) 
where M and fi are the mass and magnetic moment of 
the central object, M is the mass accretion rate and £ 
is a dimensionless constant of order 0.5-1. The funnel 
flow occurs when r m is less than the corotation radius 
r c (where the disk rotates at the same rate as the star). 
For r m <; r c , centrifugal forces may lead to ejection of 
the accreting matter ("propeller" effect). Over the last 
several decades, numerous theoretical studies have been 
devoted to understand the interaction between accretion 
disk and magnetized star (e.g., Ghosh & Lamb 1979; Aly 
1980; Lipunov & Shakura 1980; Anzer & Borner 1983 
Arons 1993; Aly & Kuijpers 1990; Spruit & Taam 1993 
Shu et al. 1994,2000; Wang 1995; Lovelace et al. 1995,1999 
Li et al. 1996; Campbell 1997; Lai 1998,1999; Terquem & 
Papaloizou 2000; Shirakawa & Lai 2002a, b; Pfeiffer & Lai 
2004; Uzdensky 2004; Matt & Pudritz 2005). There have 
also been many numerical simulations on accretion onto 
magnetic stars, with increasing sophistication and realism 
(e.g., Stone & Norman 1994; Hayashi et al. 1996; Goodson 
et al. 1997; Miller & Stone 1997; Fendt 2003; Romanova 
et al. 2003,2006; Ustyugova et al. 2006; Long et al. 2007). 

In this paper we carry out an analytical study of wave 
excitation in disks around rotating magnetic stars. The 
idea is very simple: A rotating dipole, inclined relative 
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to the stellar spin axis, induces a time-dependent mag- 
netic force on the disk. The forcing frequency normally 
equals to the spin frequency u) s , but we show that un- 
der certain conditions the magnetic force could also have 
a 2u> s component (see §2). The magnetic forces excite 
m = 1 density/bending waves in the disk, most promi- 
nently at the Lindblad/vertical resonance. This magnet- 
ically driven resonance was noted before (Lai 1999) but 
not studied. Since details of the magnetic field - disk in- 
teraction are complex (e.g., the interaction depends on 
magnetic field dissipation/reconnection), we will model 
the magnetic force using simple models that, we conjec- 
ture, capture the essential dynamical feature of the sys- 
tem. Petri (2005) examined several aspects of the forced 
oscillations in magnetized disks, but did not study the res- 
onantly excited bending waves and their evolution that are 
of interest in this paper (see also related work by Agapitou 
et al. 1997; Terquem & Papaloizou 2000). 

Wave excitations in magnetic disks may be relevant to 
understanding the flux variabilities that have been ob- 
served in many accreting systems, from pre-main-sequence 
stars (e.g., O'Sullivan er al. 2005, Bouvier et al. 2007b 
and references therein) to accreting magnetic white dwarfs 
and neutron stars. It has been suggested that the phe- 
nomenology of the qausi-periodic brightness modulations 
in CVs might be explained by magnetically excited trav- 
eling waves in the disk (Warner 2004). Of great inter- 
est are the high-frequency (kHz) quasi-periodic oscilla- 
tions (QPOs) observed in the X-ray fluxes of more than 
20 neutron stars in Low-Mass X-ray Binaries (LMXBs) 
(see van der Klis 2006 for a review). These kHz QPOs of- 
ten occur in pairs, with intriguing correlation between the 
upper- frequency v Ul the lower-frequency V\ and the spin 
frequency v s . In particular, in the accreting millisecond 
pulsar SAX J1808.4-3658, v u = 694±4 Hz, v t = 499±4 Hz, 
and v u — vi equals v s =401 Hz to within a few Hz (Wij- 
nands et al. 2003). In another accreting millisecond pulsar 
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XTE J1807. 4-294, v u — vi remains constant with an average 
value of 205 ± 6 Hz, approximately equal to v s = 191 Hz, 
even when both v\ and v u vary over a range of more than 
200 Hz (Linares et al. 2005). In some systems, however, 
v u — v\ are found to vary significantly (e.g., Circinus X-l, 
Boutloukos et al. 2006). In general, there is some evidence 
that (y u — vi)jv s clusters around 0.5 and 1 in those systems 
(about 10) for which these quantities have been measured, 
indicating that the neutron star spin plays an important 
role, although this issue is still under debate (see Yin et 
al. 2007; Mendez & Belloni 2007). It is also of interest 
to note that recent observations (Casella et al. 2007; Al- 
tamirano et al. 2007) suggest that many neutron stars in 
LMXBs have non-negligible (but complex) magnetic fields. 

Our paper is organized as follows. In §2 we discuss the 
generic forms of magnetic forces on a disk around a rotat- 
ing magnetic dipole. The basic equations governing the 
hydrodynamical response of a disk to an external peri- 
odic force are presented in §3. §4 focuses on wave dynam- 
ics around the Lindblad/ vertical resonance and corotation 
resonance and §5 gives the global solution of the disk re- 
sponse to the periodic magnetic force. Some applications 
of our result to QPOs in LMXBs are discussed in §6. 

2. PERIODIC MAGNETIC FORCES ON THE DISK 

The interaction between the dipole magnetic field of a 
star and its accretion disk is complex, and many papers 
have been published on it (§1). To focus on the dynamics 
of waves driven by magnetic forces, we consider two sim- 
ple analytical models, representing two limiting behaviors 
of the disk. The basic geometric setup is as follows: The 
disk lies in the xy-plane with its normal vector long the 
z-axis; the spin (u; s ) of the star is inclined with respect 
to the z-axis by an angle (3 and lies in the j/z-plane. Be- 
cause of the possibility of large-scale disk warping driven 
by the (static) stellar magnetic field, we will allow for 
(3 ^ 0; Lai 1999, Pfeiffer & Lai 2004; see below). The 
stellar dipole (i rotates around u> s , and the angle between 
fi and u) s is 6, such that fi varies in time according to 
fi = /i (sin 8 cos ip^x + sin 8 sin ip^y + cos 8 z) . 

2.1. Simple Model 

First consider the limiting case where the disk is a per- 
fect conductor and has no large-scale magnetic field of its 
own. The inner radius of the disk is located at r = r m . 
The magnetic field produced by the stellar dipole cannot 
penetrate the disk, and a diamagnetic surface current is 
induced. Aly (1980) has found the exact analytic solution 
for the magnetic field outside the thin disk. At a point 
(r, (f>, z = 0) (cylindrical coordinates) on the disk surface 
(z = 0, r > r m ), the field is given by 



B r 

Bio 



2fi . 



Xcos(tp- tfn) T 



Trr 3 T> 



cosx, 



— sinxsin(^- <p M ), 



(2) 
(3) 



and B z = 0, where \ — x(t) is the angle between 
fi and the z-axis. In eq. the upper (lower) sign 
applies to the upper (lower) disk surface, and T> = 

max {^Jr 2 /r^ — 1, y/2h/r n ^j , where h <C r m is the half- 
thickness of the disk. The vertical magnetic force per unit 



area on the disk is 



F z = 



2/i 2 



■ sin2x cos(ip - cp^). 
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Using the relations between various angles, we have 



(4) 



(5) 



F z =Re2_^F UJf (r)exp(i(l)-iu}ft), 

where Wf = 0, ±w s , ±2lu s , with 

F (r) =-iF D {2 -3 sin 2 8) sin 2(3, 
F Us (r) = T D sin 28 (cos (3 + cos 2/3), 
F_ Ws (r) = T D sin 28 (cos (3 - cos 2/3), 
Fiu> 3 ( r ) — sin 2 8 sin (3 (cos (3+1), 

F-iw a ( r ) — sin 2 8 sin (3 (cos [3 — 1), (6) 



and 



•F D 



(7) 
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The disk is unlikey to be completely diamagnetic: var- 
ious instabilities and dissipation mechanisms will allow 
some of the stellar magnetic field to partially penetrate 
the disk. As a simple model, we decompose the vertical 
stellar dipole field on the disk into a static component and 
a time- varying component: 

= — — cos x = (sin (3 sin 8 sin uj s t — cos (3 cos 8) . 

(8) 

We assume that the static field penetrates the disk, while 
the the variable field is shielded out of the disk by a screen- 
ing current. Because of the shear between the disk and 
the plasma outside the disk, the threaded vertical field is 
wound to produce an azimuthal field that has different 
signs on the upper and lower surfaces of the disk. We thus 
adopt the following ansatz for the magnetic field in the 
disk (see Lai 1999) 

2[i 4/i 

— sinxcos(a3 — Wu ± — — sin (3 sin 8 sin u> s t, (9) 
r 6 txt^u 

B v = -^■sinxsin(( y 9 — iy9 M )±C-^-cos/3cos0, (10) 
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cos (3 cos 8, 
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where £ is a positive constant of order unity. The vertical 
magnetic force (per unit area) on the disk is given by 

4^ 2 

F z = sin/3sin#sinw£sinxcos(i/2 — ip^) 
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cos (3 cos 8 sin x sin(ip — ip^). (12) 



This force can be written in the form of eq. ([5]), with 
F (r) = iF D sin 2 8 sin 2(3 + F T cos 2 8 sin 2(3, 



F Ue (r)=J : D sin 28 sin 2 (3 + iT T sin 28 cos (3 cos 2 | , 

F-uj 3 (r) = —J~d sm 20 sin 2 (3 + iTt sin 28 cos (3 sin 2 — , 

F 2Us (r) = -iT D sin 2 8 sin (3 (2 + cos (3) , 
F-2lo s (r) = iJ : D sin 2 8 sin (3 (2 — cos /3), 
where 

Cm 2 



4irr 6 ' 
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2.2. General Property of Magnetic Forces 

The magnetic forces discussed above share some general 
properties. The zero-frequency force components F have 
been shown to lead to secular warping and precession of 
the disk (Lai 1999; Shirakawa & Lai 2002a,b; Pfeiffer & Lai 
2004): The force resulting from threaded field (oc Tt) in- 
duces a warping instability, while the force resulting from 
the dielectric response of the disk (oc Tu) gives rise to pre- 
cession. The negative-frequency components {F— Ug and 
F~2u s ) average to zero on the timescale of spin period. 

The positive- frequency forces, F UJb and F% Ull , are of in- 
terest for this paper. They are generally present for non- 
zero magnetic field inclination angle 8. Note that the 2oj s 
component arises because the dipole field varies as sinw s f 
or cos u) s t, and the induced screening current also varies as 
sinw s t or coso; s i, and F^, is present only when sin/3 ^ 0. 

We note that periodic magnetic forces in the radial di- 
rection are also present in some disk models, but we we 
will focus on the dynamics of the disk under the periodic 
vertical force in this paper. 

3. BASIC EQUATIONS: DISK DYNAMICS UNDER A 
PERIODIC NON-POTENTIAL FORCING 

The key approximation underlying our study is that we 
treat the waves in the disk using pure hydrodynamics, with 
the magnetic effect entering only as external forces. This 
simplification allows us to focus on the dynamics of wave 
excitation, propagation and dissipation. As we will see, 
the wave dynamics is sufficiently complex and subtle even 
without magnetic fields, and we believe that such a hy- 
drodynamical treatment can be considered a useful first 
step. Future study using MHD will be of great interest 
(and obviously much more difficult). We conjecture that 
in a real system, the wave dynamics will not differ qual- 
itatively from that described in this paper since the the 
waves of interest here are basically sound waves modified 
by disk rotation, and in the presence of magnetic fields 
the sound waves simply become fast magnetosonic waves 
(with the sound speed replaced by fast wave speed). 

In this section we develop the general hydrodynamical 
equations for three-dimensional (3D) disks under generic 
non-potential forces (cf. Zhang & Lai 2006). We con- 
sider a geometrically thin gas disk and adopt cylindrical 
coordinates (r, ip, z). The unperturbed disk has velocity 
vo = (0,rf2, 0), where the angular velocity f2 = f2(r) is 
taken to be a function of r alone. The disk is assumed 
to be isothermal in the vertical direction and non-self- 
gravitating. Thus the vertical density profile is given by 



Po{r,z) 



exp(-Z 2 /2), with Z = z/h (15) 



where h — h(r) — c/Cl± is the disk scale height, c = c(r) = 

y/po/po is th e isothermal sound speed, a = cr{r) — J dz po 
is the surface density, and Oj_ is the vertical oscillation 
frequency of the disk. 

We now consider perturbation of the disk driven by an 
external force (per unit mass) f . For simplicity, we shall 
assume that the perturbation is isothermal, so that the 
(Eulerian) density and pressure perturbations are related 
by SP = c 2 Sp. The linear perturbation equations read 

^ + (v • V)u + (u • V)v - -Vr) + f , (16) 



dp 

-£ + V • (p u + v 5p) = 0, 
at 



(17) 



where u = 5w is the (Eulerian) density perturbation, and 
rj = dp/po is the enthalpy perturbation. Without loss 
of generality, each perturbation variable and the external 
force are assumed to depend on t and ip as e~xjp(imip — iuit), 
where to > is an integer and u) is allowed to be either 
positive or negative, corresponding to the prograde or ret- 
rograde wave, respectively. Equations (fTBj) and (TTT]) then 
reduce to 

d 



2Qu u 



c). 



-V + fr, 



k im 
iuu v + —u r = — —77 + f v , 



2n 



(18) 
(19) 
(20) 
0(21) 

(22) 
(23) 

In this paper we will consider cold, (Newtonian) Keplcrian 
disks, for which the three characteristic frequencies, il, f2j_ 
and k, are identical and equal to the Keplerian frequency 
Vl K = (GM/r 3 ) 1 / 2 . However, we continue to use differ- 
ent notations (f2, Q j_, k) for them in our treatment below 
when possible so that the physical origins of various terms 
are clear. 

To separate out the z-dependence, we expand the per- 
turbations with Hermite polynomials H n (cf. Okazaki et 
al. 1987; Kato 2001) 

' frn(r) ' 



d 

-IUJU Z = --q^V + Jz, 

_ Po 1 9 , - im d 

-iui^rr]+ - T rp !i r -I p u v + -T-(p u z 

& r or r oz 

Here ui is the "Doppler-shifted" frequency 

lu = lu — mCl, 

and k is the epicyclic frequency, defined by 

K 2 = ^4(r 2 tt). 
r dr 



fr(r,z)' 
U(r,z) 
■q(r, z) 
u r (r,z) 
u v (r, z)_ 

fz{r,z) 
u z (r,z) 



E 



E 



ftpn 

Vn(r) 
u rn (r) 

Uir) 
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(24) 



where H' n — dH n /dZ, and the Hermite polynomial is de- 
fined by H n (Z) = (-l) n e z2 / 2 d n (e~ z2 / 2 )/dZ n . Note that 
since Hi = Z, H2 = Z 2 — 1, the n = 1 mode coincides with 
the bending mode studied by Papaloizou & Lin (1995) 
(who considered disks with no resonance), and the n = 2 
mode is similar to the mode studied by Lubow (1981). 
With the expansion in ([24]) . the fluid equations (|2Tj) be- 
come 

np 



-iCbUm - 2flu vn = —-j-r]n H V 
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(26) 
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where /i = d In h/d In r. Eliminating u vn and u zn from 
eqs. (|25 ]) -([28 |) . we have 
efyn 2mf7 D . a 

-j— = — —Vn ~ —iu rn + - w/„ + (n + l)(n + 2)r] n+2 \ 
dr rw uj r 



.20 

~frn ~t~ ^ ~ fipm 



(29) 



c£it r 



<iln(r<7) tok 2 



2rO£ 



h 2 



rVn (nu rn + u r ,„_ 2 ) + — Un + -rr-fzn, (30) 



cr r 
where we have defined 

D = k 2 - Cj 2 



K 2 -(w-mS]) 2 . (31) 
We can check that for an external potential force, f = 
— V(/>, these equations reduce to those of Zhang & Lai 
(2006). 

Consider local free wave solution of the form rj n oc 
exp[i fk{s)ds]. For \kr\ > 1, from the eqs. ( ^ -([gD ]) . 
in the absence of the external force, we find (see Okazaki 
et al. 1987; Kato 2001) 

(w a -K a )(w a -nfi x )/£ a = fc a c 2 , (32) 
where we have used h = c/f2j_ <C r (thin disk), and 
m, ?i -C r//j. The dispersion relation is useful for iden- 
tifying wave propagation zones in the disk (see Zhang & 
Lai 2006). 

4. WAVES AT LINDBLAD/VERTICAL RESONANCE AND 
COROTATION RESONANCE 

We now consider the magnetic force on the disk as given 
in §2, i.e., the force per unit area is F z — F u (r) exp(ic/j — 
iwt), with uj — uj s or 2lu s . This corresponds to m = n = 1, 
with 

fzi(r) = -F u (r), f r = f =0, (cj = uj s ,2uj s ). 
a 

The wave propagation diagram [based on the dispersion 
relation (|32p with n = m = 1] is shown in Fig. 1. 

Note that for w < = k = f2xj the m = n = 1 mode 
dispersion relation (|3"2"|) reduces to u = ±cfc/2, i.e., the 
bending wave propagates non-dispersively at speed c/2. 
This is the regime explored in a number of previous works 
(e.g., Papaloizoufe Lin 1995; Ogilvie 2006). However, here 
we are particularly interested in the disk region where f2 
is of the same order as uj. 




Fig. 1.— A sketch of the function G = (k 2 - tj 2 )(l - f^/iS 2 ) 
as a function of r for m = n = 1. The WKB dispersion relation is 
G = —k 2 c 2 , and thus waves propagate in the region G < 0. Note 
that the Lindblad/ vertical resonance (L/VR) is the turning point 
k — > and the corotation resonance (CR) is a singularity (k —* oo). 



Consider eqs. (l2"9t (|3T))) for n = 1. The coupling terms 
can be neglected if ^ \rfa\, which is justified since there 
is no n = 3 driving force. Keeping only the potentially sin- 
gular terms (oc D~ 1 , Gj~ 2 ), we have 
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(34) 



This equation can be simplified further by changing vari- 
able via 

m (r) = (-D/ra) 1/2 y(r). (35) 
Keeping the leading-order terms, we obtain 

' d 2 (uj 2 - n 2 ± )D 3(Cjdn/dr + ndn/dr) 2 
dr 2 



h 2 uj 2 nl 

(-roD) 1 ' 2 



D' 2 



huj 



7,2 



fz 



(36) 



The terms dropped inside [...] in the above equation are 
considered to be smaller than the retained terms for the 
following reasons: Firstly, away from the D = region, 
they are negligible compared to the term cx h~ 2 owing 
to the assumed smallness of h, and secondly, they are in- 
significant with respect to the term cx D~ 2 in the vicinity 
of D = 0. 

Equation (|36p reveals two special points: one is at uj = 



(corotation resonance), the other is at uj 



(Lind- 



blad resonance). The corotation resonance is a singu- 
larity of the wave equation, and the singularity implies 
that a steady emission or absorption of waves may oc- 
cur there. The Lindblad resonances are turning points 
at which waves are generated or reflected. Another in- 
teresting point is at uj 2 = fi^. It is a transition point 
where the vertical resonance occurs. For a Keplerian disk, 
k = fl± = £1, the vertical resonance coincides with the 
Lindblad resonance, and we will call it Lindblad/ vertical 
resonance (L/VR). Equation (f3"6"|) then reduces to 

~ <P_ ^ (uj 2 — Vi 2 ) 2 _ 3uj 2 (dn/dr) 2 ' 
'dr 2 ^ 



h 2 



-Cj 2 n 

a^r? (uj 
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i/~2 



(UJ 2 - K 2 ) 2 
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huj 2 

This is our basic working equation. 



fzi = S u {r). (37) 



4.1. Lindblad /Vertical Resonance 

We now study wave excitations near a Lindblad/ vertical 
resonance (L/VR), where D = 0, or uj = k = il, uj — 20, 
and the (outer) L/VR radius is denoted by rr,. Changing 
r to the new variable x = (r — Tl)/tl an d keeping the 
leading order terms, we find that, for \x\ <ti 1, eq. ()37D 
reduces to 

( b 2 x 2 -—Av = S u rl, (38) 



dx 2 
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Ax 2 



where 



b 2 = 



64 (dVt/dr) 2 r A 
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= 36 
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(The subscript "L" implies that the quantity should be 
evaluated at r = Vl). The two independent solutions of 
the homogeneous version of eq. psp are 

y± cx (bx)-V 2 e ±ibx */ 2 . (40) 
The general solution of the inhomogeneous eq. (f3"5| is then 



y 



-y-\ 



V-S u rl f x y+S„rl „ 

■ ax+y- J — — — ax+C+y^ 



W 



w 



-C-V-, 
(41) 

where W = y+dy- / dx — y_dy+/dx = —2i is the Wron- 
skian. The constants C± can be fixed by requiring waves 
propagating away from the resonance. Using eq. (|35p. we 
obtain 



rh \ 1 ^ 2 



< 2 /2 



c 



dc 



-zC 2 /2 



dc 



(42) 



where / zl is evaluated at r = tl, and we have defined 
C = b x / 2 x. Note that although our analysis here is lim- 
ited to \x\ -C 1, we have extended the integration limit to 
C = ±oo in eq. (|42|) . This is valid because b = Gr^/h ^> 1 
for a thin disk and the integrands in the integrals are highly 
oscillatory for |£| ^> 1 (so that the contribution to the in- 
tegrals from the |£| 3> 1 region is negligible). 

To calculate the angular momentum transfer through 
the resonance, we note that for £ — > +oo (but still i<l), 
equation (|42"]) becomes 

The angular momentum flux to the r > Tl is 



oo). (43) 



lira ( drjT 
F (r>r L ) = —Im\r ll — 



-\fz 



^(rVl/.l 



2 \dD/dr 
Similarly, for £ — > — oo, equation (142|) gives 



.(44) 
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1/2 



= i - /, 1 (2«) 1 /2 e -C/2 („--«,). (45) 



The angular momentum flux to the r < rt region is iden- 
tical to F(r > r£). The total torque on the disk acted 
through the OL/VR is then 

T L = 2F(r>r L ) = ^{r 2 a\f zl \ 2 ) L . (46) 

To estimate the torque Tl, we consider u> — u> s , and use 
fzi = F^/cr with F u ~ ^ 2 /(47rr 6 ) (see §2). For a Kep- 
lerian disk, Vl = 2 2 / 3 r c and we assume r c /r m = A (for 
accretion to occur, we require A > 1). Using eq. (1), we 
then obtain 

M 2 

r £ ~4xnr 4 (47) 

(£A) 7 cr(r L ) 

The canonical accretion torque is Ta = M{GMr m ) 1 / 2 . 
Using M = 2tt (ra\u r \)L (where u r is the radial velocity if 
the unperturbed accretion flow), we find 

(48) 

Thus the resonant torque is much smaller than the accre- 
tion torque. 



^~3x 1(T 3 £- 7 A- 6 - 5 
Ta 



4.2. Corotation Resonance 

The corotation resonance (CR) is located where Co = 
or lo = mVt{r c ). The WKB dispersion relation shows that 
for n = 0, waves are evanescent in the region around the 
corotation radius r c , while for n > wave propagation is 
possible around r c . 

We again focus on the n = m = 1 waves. In the vicin- 
ity of corotation, the terms oc h~ 2 in eq. (f3~4"|) or (j3"T|) are 
dominant, and we only need to keep these terms and the 
second-order differential term. Equation (j3~4"|) then reduces 
to 

rl 2 D{Q 2 -n 2 ) D 

-Vi = —hi- ( 49 ) 



dr 2 



h 2 Co 2 nl 



Defining x — (r 
x = 0, we have 
d 2 



hCo 2 ' 

r c )/r c and expanding eq. (|49|) around 



C 



dx 2 (x + ie) 2 



= C 



hf z 



where hf z \ is evaluated at r = r. 
C = 



(*- 
and 

2r x 2 
3^ 



ie) 2 



» 1 



(50) 



(51) 



/i e?£7 1 dr t 

(where the subscript "c" means that the quantity is evalu- 
ated at r = r c ). In eq. (|50p . we have inserted a small imag- 
inary part it (with e > 0) in l/x 2 because we consider the 
response of the disk to a slowly increasing perturbation. 
The general solution to eq. (|50| is 

Vi = hf z i + Mz x l 2 z iv + Nz x l 2 z- lv 

= hf zl + Mz 1 ' 2 e ivlnz + Nz 1/2 e- iulnz , (52) 



where 



i/=yC-|»l, z = x + ie (with e > 0) 

and M and are constants. The first term, the non- 
wave part, is a particular solution, while the other two 
terms are solutions to the homogeneous equation, depict- 
ing the waves. The z 1 ^ 2 z lu term has a local wave num- 
ber k = d(vhiz) I dr = v/(r c x), with the group velocity 
v g = dtu/dk — —Q/k — —3r c x 2 Q c /(2v) < 0, thus it rep- 
resents waves propagating toward small r. Similarly, the 
z v l 2 z lv term has v g > and represents waves propagating 
toward large r. 

As shown in Zhang & Lai (2006), waves with n = 1 
can propagate into the corotation region and get absorbed 
there. Consider an incident wave propagating from the 
x > region toward x = 0: 

m {x > 0) = A+x 1/2 e™ lnx , (53) 

where A + is constant specifying the wave amplitude. The 
transmitted wave is given by 

rji {x < 0) = i A+ e-™ {-xf/ 2 e iv ln( ~ x) . (54) 

Since v 3> 1, the wave amplitude is vastly decreased by a 
factor e~^ v after propagating through the corotation. The 
net angular momentum flux absorbed at corotation is 

where in the last equality we have used v ~ \/~C 3> 1. 
Thus, a wave propagating from r > r c into the corotation 
resonance deposits almost all of its positive angular mo- 
mentum at r = r c . Similarly, a wave propagating from 
r < r c into the corotation resonance deposits its negative 
angular momentum at r = r c . 
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5. GLOBAL DISK RESPONSE 

Having studied the behavior of waves near L/VR and 
CR, we can now construct global solution for the waves 
excited by the external force. 

Away from the L/VR, the two linearly independent 
WKB solutions to the homogeneous wave eq. (|3T[) are 



k r dr 



where 



V± = K 1/2 exp (± i J 

- 9 - ft 2 _ uj{uj - 2ft) 



/iftcj 



hVL(uj - ft)' 



(56) 



(57) 



It is easy to check that eq. (I5"I51) reduces to (|4T))) for 
\r — tl\ thus it is valid even near the L/VR. From 

cq. (j3"5j) , we have 

'/ift(w-ft) ll/2 



Vi± = 



TO 



exp ± i / k r dr 



(58) 



Near the CR, this solution reduces to 
m± oc x 1/2 e Twlnx , for |x| = |(r - r c )/r c \ < 1. (59) 

Thus equation (|58|) represents the two independent solu- 
tions of the homogeneous wave equation for all radii. 

As discussed before, Waves are mainly excited at the 
L/VR, where waveform is given by eq. (|4"/T)) . Matching 
eq. (|4*3")l or (|4"5|) with the general solution, r]i = C + t]i + + 
C-r\\- (where C± are constants), we find that the waves 
away from the L/VR are: 



m = 



■ 2 \ 1/2 



3lo 



L 

x exp | / 



■ 2 \ I/ 2 



fzi{rh) 
k r dr 



Ml(u - ft) 



1/2 



3w 2 



x exp 



fzi(r L ) 



(r > r L ) 
/ift^-ft)" 1172 



(60) 



i J k r dr 



rcr 

(r c <r < r L ).(61) 



Note that although for r c < r < Tl the wave (|6T|) has 
positive phase velocity (since k T < for r c < r < r^), the 
group velocity is negative, i.e. the wave propagates toward 
small radii. As this wave approaches r c , the waveform re- 
duces to [see eq. ([55])] 



% = A+x 1 ' 2 e wlnx , 
where v = 2r c / (3h c ) and 

A+ = 



0<x = - — - < 1 



(62) 



2 ^T 



l"l"rr) ; f z i(r L ) 



- 

2ra 



x exp 



k r dr — iv In x+ 



(63) 



where < x + = (r + — r c )/r c <C 1. After passing through 
the CR, the wave amplitude is significantly reduced: 

iA + e-™(-x) 1/2 e ivln( -- x \ 
r - r r . 



m 



<0, x«l 



(64) 



This then joins onto the inward-going wave solution 



m 



3uj 2 



1/2 



fz\(r L ) 



htt{uj - ft) 



1/2 



x exp 



k r dr 



(r < r c ), (65) 



with <p — J r + k r dr + i/ln(— x+/x_), where x_ = (r_ — 
Tc)/r c < and |x_| -C 1. Using eq. (|55|) . we find that the 
angular momentum flux deposited at r c is 

AF r = ^(r 2 a\f z i\ 2 ) T , (66) 



3w 2 v "U^l /£> 

in agreement with the angular momentum flux emitted 
from toward corotation [see eq. (|46|) ]. 

Figure 2 depicts the global solution describing waves ex- 
cited at the L/VR, which either propagates outwards or 
inwards toward the CR at which the angular momentum 
deposition occurs. 




Fig. 2. — The enthalpy perturbation rji of the wave driven by 
an external force (m = n = 1). The wave is excited at the L/VR 
(r/r^ = 1), and propagates toward large radii and small radii. The 
inward-going wave is absorbed at the corotation resonance. The 
lower panel focuses on the region between corotation and rj, . 



6. APPLICATION TO QPOS 

The main result of the paper is illustrated in Fig. 2. 
It shows that under a periodic magnetic vertical forcing, 
the disk responds by launching m = n = 1 waves from the 
Lindblad/vertical resonance (where u> = 2ft, with the forc- 
ing frequency uj = oj s or 2w s ). The wave either propagates 
to large radii or propagates inward where it gets absorbed 
at corotation resonance (where u> = ft). 

We have already shown that the torque carried by the 
excited wave is small compared to the canonical accre- 
tion torque (see §4.1). Nevertheless the wave may mani- 
fest itself by inducing or modulating variabilities. Clearly, 
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the wave is most visible at the Lindblad/ vertical reso- 
nance since the wavelength is the largest there. Vari- 
ations caused by the wave are most likely due to fluid 
elements around rj,. From eqs. (|oT))) - (|61[) . wc find that 
the amplitude of enthalpy perturbation at rj, is given by 
1 771 1 = (■7Tr/12h) 1 ^ 2 (F UJ h/a) (all quantities are evaluated at 
r = r L ). With F u — /a 2 /(4-7rr 6 ) (see §2), r m = r c /X (where 
A > 1 is a constant), — 2 2 / 3 r c , and using eq. (1) and 
M = 2nr<j\u r \, we find the dimensionless amplitude 

^^0.05 ( AO- 7 / 2 (g)^, (67) 

where u r is the radial (inflow) velocity. For a a-disk, 
\u r \ ~ (3/2r)ahc not too close to the disk inner edge, 
we have 

M. .f«(Aer 7 / 2 (^ /2 . (68) 

The dimensionless vertical velocity perturbation, |u zl |/c, 
is the of the same order of magnitude at r ~ tl as 1 771 1 /c 2 . 
Since A£ ~ 1, the dimensionless perturbation amplitude 
may reach a few percent, while definitely remaining in 
the linear regime. Note that the vertical velocity per- 
turbation amplitude increases toward the corotation ra- 
dius, non-linear saturation will occur there. But since the 
wavelength approaches zero as r — ► r c , variations caused 
by fluid elements around corotation resonance may not be 
visible. 

It seems inevitable that bending waves studied in this 
paper will be excited in disks around magnetic stars. We 
may suggest various ways that these resonantly excited 
waves may induce variabilities and QPOs. In the context 



of LMXBs (van der Klis 2006), one may imagine that the 
oscillating fluid perturbation at the Lindblad resonance 
can produce a beat phenomenon by modulating the "seed" 
radiation from the inner region of the disk (e.g. Lamb & 
Miller 2003). Suppose the "seed" radiation has a quasi- 
periodic variability with frequency vh (e.g. due to orbital 
motion of blobs at the disk inner edge or some other mech- 
anisms). The fluid element at the Lindblad/ vertical reso- 
nance has orbital frequency il = u>/2 (since lu — Q = Q± 
for m = 1, and fl± ~ f2). The reprocessed radiation would 
then show an additional QPO at frequency v\ = Vh — u/2. 
In this regard, it is particularly interesting that in our 
model, the driving frequency v — u>/(2tt) can be either 
the spin frequency v s or twice of that, depending on the 
disk geometry with respect to the stellar spin and dipole 
axes (see §2). As noted in §1, both Vh — vi = v s and 
Vh — v\ = v s j2 have been observed in LMXBs. This fea- 
ture has not been explained by current models (see van 
der Klis 2006). For many systems, the beat is not per- 
fect. This might be due to the fact the excited wave still 
has significant strength away from the Lindblad/ vertical 
resonance (see Fig. 2). 

Finally we note that our treatmnet of wave excitation 
by magnetic forces is performed under simplifying assump- 
tions, as detailed in section 3. We hope that the novel fea- 
tures of our model in connection with QPOs will motivate 
future, more rigorous studies of the dynamics of waves in 
disks around magnetic stars. 



This work has been supported in part by NASA Grant 
NNX07AG81G and by NSF grant AST 0707628. 
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